UNCLASSIFIED 


Defense  Technical  Information  Center 
Compilation  Part  Notice 

ADPO 11592 

TITLE:  Homogenization  Formalisms  for  Nonlinear  Particulate  Composite 
Mediums 

DISTRIBUTION:  Approved  for  public  release,  distribution  unlimited 


This  paper  is  part  of  the  following  report: 

TITLE:  International  Conference  on  Electromagnetics  of  Complex  Media 
[8th],  Held  in  Lisbon,  Portugal  on  27-29  September  2000.  Bianisotropics 
2000 

To  order  the  complete  compilation  report,  use:  ADA398724 

The  component  part  is  provided  here  to  allow  users  access  to  individually  authored  sections 
of  proceedings,  annals,  symposia,  etc.  However,  the  component  should  be  considered  within 
the  context  of  the  overall  compilation  report  and  not  as  a stand-alone  technical  report. 

The  following  component  part  numbers  comprise  the  compilation  report: 

ADP01 1588  thru  ADP01 1680 


UNCLASSIFIED 


23 


Homogenization  Formalisms  for  Nonlinear 
Particulate  Composite  Mediums 

A.  Lakhtakia1.  M.  N.  Lakhtakia2  and  W.  S.  Weiglhofer3 


1 CATMAS  — Computational  & Theoretical  Materials  Sciences  Group 
Department  of  Engineering  Science  and  Mechanics,  Pennsylvania  State  University 
University  Park,  PA  16802-6812,  USA 

Tel:  +1  814  863  4319;  Fax:  +1  814  863  7967;  email:  AXL4@psu.edu 

2 1811  Red  Lion  Drive,  State  College,  PA  16801-3012,  USA 

3 Department  of  Mathematics,  University  of  Glasgow 
University  Gardens,  Glasgow  G12  8QW,  Great  Britain 

Abstract 

The  Maxwell  Garnett  (MG)  and  the  Bruggeman  (Br)  formalisms  are  extended  to  homogenize 
nonlinear,  two-component,  composite  mediums.  The  chosen  material  topology  is  ellipsoidal 
and  weak  nonlinearity  is  assumed.  The  MG  formalism  is  illustrated  by  a case  in  which 
both  component  materials  are  bianisotropic,  but  only  the  inclusion  component  is  nonlinear. 

Both  component  materials  are  isotropic  dielectric  in  the  following  case,  but  only  one  has 
an  intensity-dependent  permittivity  scalar;  and  the  Br  formalism  is  applied  to  show  that 
the  homogenized  composite  medium  is  anisotropic  and  has  cubically  nonlinear  dielectric 
properties.  Enhancement  of  nonlinearity  emerges  as  a significant  possibility. 

1.  Introduction 

The  Maxwell  Garnett  (MG)  and  the  Bruggeman  (Br)  formalisms  for  the  homogenization  of  linear 
composite  mediums  formed  by  randomly  dispersing  electrically  small,  bianisotropic,  ellipsoidal 
inclusions  in  bianisotropic  host  materials  can  be  extended  to  composite  mediums  comprising 
nonlinear  materials.  We  illustrate  this  extension  with  an  example  for  each  formalism.  The 
nonlinearity  is  assumed  as  weak  and,  therefore,  perturbatively  tractable. 

In  the  first  example,  both  component  materials  are  bianisotropic,  but  only  the  inclusion  ma- 
terial is  nonlinear.  The  effective  constitutive  properties  of  the  homogenized  composite  medium 
(HCM)  are  estimated  using  the  MG  formalism.  The  linear  and  nonlinear  properties  of  the  HCM 
are  estimated  separately  in  consequence  of  two  assumptions:  the  nonlinearity  of  the  inclusion 
material  is  weak,  and  the  composite  medium  is  dilute  [1]. 

In  the  second  example,  both  component  materials  are  isotropic  dielectric,  but  just  one  of 
them  has  an  intensity-dependent  permittivity  scalar.  Application  of  the  Br  formalism  shows 
that  the  HCM  is  anisotropic  and  has  cubically  nonlinear  dielectric  properties.  The  anisotropy  of 
nonlinearity  can  be  considerably  different  from  the  anisotropy  of  the  linearity;  and  the  possibility 
of  nonlinearity  enhancement  exists  [2], 

2.  Nonlinear  Bianisotropic  Composite  Medium 

Consider  a countably  infinite  number  of  identical,  electrically  small,  ellipsoidal  inclusions  that 
are  similarly  oriented  but  randomly  embedded  in  a host  material.  Each  inclusion  has  a volume 
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v,  the  number  density  of  inclusions  is  denoted  by  JV;  while  / = Nv,  0 < / < 1,  is  the  volumetric 
proportion  of  the  inclusion  medium.  Typically,  / < 0.2  for  the  MG  formalism  to  yield  adequate 
results.  The  inclusions  are  described  by  a shape  dyadic  which  is  real  symmetric  with  positive 
eigenvalues  [3]. 

The  use  of  6-vectors  and  6x6  dyadics  is  very  convenient  for  bianisotropic  mediums  as  it 
permits  a compact  notation.  In  this  notation,  the  constitutive  properties  of  the  host  material 
are  expressed  as 

®“> = ( ft!  ft) ) ■ EH = ■ EM  ■ (l) 

where  the  6-vectors  F(w)  = [E(uj),  H_{ u)]T  and  G(uj)  = [I)(u;),  B(u)]T,  the  superscript  T 
denoting  the  transpose. 

With  the  6x6  constitutive  dyadic  of  free  space  denoted  by  Qv  (which  contains  eo  and  /io), 
the  linear  and  the  nonlinear  constitutive  properties  of  the  inclusion  material  are  best  expressed 
through  the  6-vector  Q(u>)  = [P{u>) , M(u))]T  = G(u;)  — • F(a;),  which  contains  both  the 

polarization  field  P{w)  and  the  magnetization  field  M(lo).  This  6- vector  is  split  into  linear  and 
nonlinear  parts  as  Q(cj)  = Q*(u;)  + Qnl(u).  Its  linear  part  Q?(oj)  obeys  the  constitutive  relation 

Qf(cd)  = Cm(o;)  — Qv  • F(w),  with  Qtn(cj)  analogous  to  Qk(uj). 

The  nonlinear  properties  of  the  inclusion  material  are  described  as  follows:  Let  W = 
w3, . . . ,u>m}  be  a set  of  M > 1 angular  frequencies,  there  being  no  requirement  that  all 
members  of  W be  distinct.  Suppose  there  exists  an  ensemble  of  M fields  F(u;m),  (1  < m < M ). 
Then  the  simultaneous  action  of  this  ensemble  of  fields  in  the  inclusion  material  gives  rise  to  the 
nonlinear  part  of  Q(u»)  at  u)  — u )nl.  The  j- th  element  of  this  6- vector  is  given  by 

6 6 6 6 f M 

o?V‘)=  E E - E - E 

jl=l  J2  = l jm  — 1 JM=1  I "=1 

where  (^ni;  W)  represents  the  nonlinear  susceptibilities  of  the  inclusion  material. 

The  angular  frequency  unl  is  simply  related  to  all  members  of  W as  ujnl  = oi  uji  + 02  W2  + • • ■ + 
bmwm,  with  am  = ±1,  (1  < rn  < M)\  furthermore,  u>nl  may  or  or  may  not  lie  in  W.  When 
an  = —1,  Fjn( ajn)  must  be  replaced  by  its  complex  conjugate  on  the  right  side  of  (2)  and  in 
subsequent  derivative  expressions. 

Details  of  the  implementation  of  the  MG  formalism  to  homogenize  the  described  composite 
medium  were  given  elsewhere  [1].  In  summary,  the  constitutive  relation  of  the  HCM  is  given  by 

G(u;)  = Qmg(oj)  • E(oj)  + s(u, u>nl)  Qn*s(u>) , u)  € (w  U {a;nZ})  , (3) 

where  the  switching  function  s(ui,u}nl)  equals  unity  when  both  of  its  arguments  are  the  same, 
but  is  null-valued  otherwise.  Expressions  for  Cmg(oj)  have  been  available  for  about  three  years 
[3],  and  need  not  be  reproduced  here.  The  nonlinear  source  polarization-magnetization  field  is 
expressed  through  [1] 

Qnls(unl)  = 1 1 - iunl  l - iojnlfUm/h{ojnl)  • ain/h(u)nl)  1 • Din/ft(a;nJ)|  • Qn*ff(unl) , (4) 
where 

(«}/)  (u."1)  = / EE-  E - E 

il  = lj2=l  jm=  1 jM  = 1 I 

Mr  'J 

1 1 [(xin/MG(^n)  • EK)).J  1 , 1 < j < 6 . 


j,  1 < i < 6 , 


(5) 
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In  these  expressions,  the  6x6  dyadics 

iin//l(o ;)  = [gin(u>)  - Qh(u})\  • Xin/h(w) , (6) 

Xin/PM  = {i  + iujUin/pH  • [giT»  - SPM]  }_1  , p = h,  MG,  (7) 

while  Dm/p  (resp.  Dm/,p)  is  the  6x6  depolarization  dyadic  of  an  ellipsoidal  (resp.  spherical) 
exclusion  region  in  a linear  medium  with  Qf  as  its  constitutive  dyadic  [3]. 

The  presented  formalism  is  general  in  that  it  can  be  used  to  examine  harmonic  generation, 
parametric  oscillation,  self-focusing,  stimulated  Raman  scattering,  and  a multitude  of  nonlinear 
phenomenons.  It  also  updates  and  extends  our  previous  work  on  the  MG  formalism  for  complex 
nonlinear  composite  mediums  [4,  5]. 


3.  Anisotropic  Dielectric  Composite  Medium  with  Intensity— Dependent  Per- 
mittivity Dyadic 

Next,  we  implement  the  Br  formulation  for  homogenizing  a mixture  of  two  dielectric  materials. 
For  the  sake  of  illustration,  here  only  one  component  material  is  assumed  to  be  nonlinear:  it 
possesses  an  intensity-dependent  permittivity  scalar.  Accordingly,  gh  = eh  /,  gin  — ein[_  = 
(e£n  + e|J}|i£|2)  I,  fih  = ptn  = /iol,  and  the  remaining  components  of  Qh  and  are  null- 
valued. The  ^-dependences  are  not  explicitly  identified,  as  both  the  linear  and  the  nonlinear 
fields  vibrate  at  the  same  frequency.  Both  component  materials  are  assumed  to  have  parallel 
ellipsoidal  topologies.  As  both  component  materials  are  treated  in  the  same  manner  in  the  Br 
formalism,  the  labels  in  and  h lose  the  meanings  they  have  for  the  MG  formalism,  and  the  results 
are  valid  prima  facie  for  / € [0,  1], 

The  Br  formalism  requires  the  solution  of  the  dyadic  equation  [3] 

/ (ein£-  gBr)  • (iin/Br)_1  + (l-/)  (eftI-£Br)  * (&h/Br)  (8) 


where 


Xfn/Br  gin/Br  . ^ in  £ _ ^BrJ 


(9) 


and  is  defined  similarly;  while  gn/Br  and  ^/Br  are  3x3  depolarization  dyadics. 

The  HCM  is  anisotropic,  and  the  Br  formalism  predicts  its  permittivity  dyadic  as  [3,  2] 


gBr  = eBr  + eBr  uyuy  + eBr  uzuz . 


(10) 


A perturbative  treatment  permits  the  ansatz  gEr  ~ g£r  + e^r  E , consistently  with  our  as- 
sumption that  the  nonlinearity  of  e,n  is  weak;  hence, 


eBr  ~ eBr  -j-  eJ 

— CX<  — cXn/ 


.Br 


E 


. Br  ~ Br  , Br 
y ~ cyi  ~ cj/ni 


E 


.Br  ~ Br  , Br 

* — Zzi  ~ CZ„( 


E 


Therefore,  the  Taylor  expansions 

£p/Br  ~^/Br+£p/Br'  12 


E 


p — in,  h. 


(11) 


(12) 


emerge.  Expressions  for  P^Br  are  available  in  Ref.  [2].  Accordingly, 


£p/Br~£p/Br+^{Br 


E 


Br 


p = in,  h, 


(13) 


26 


where 


gn/Sr  = j + ^ gn/Br  . | _ grj,  (14) 

^/Bt  = J_  + iuj  £j/Br  • (ehL-  |f r)  , (15) 

l£/Br  = » [£T/Br  • (» 43  I - (5r)  + a^/jBr  • (4"  l - *r)] ' (>«) 

^Br  = w [-  a?/Br  • sf,r + • («*  i - if')] . (i7) 

The  local  field  factor  g is  estimated,  as  a first  approximation,  as  g (l/9)|trace  |^n^Br  J |2 
where  Z?/Br  = (^n/Br)_1  [2], 

With  the  foregoing  developments,  the  nonlinear  dyadic  equation  (8)  separates  into  two  parts 

[2]:  (i)  |fr  is  the  solution  of 

ifr  = [/  x?/Br  + (i  - /)  £n/Br] _1  • [/4n  2$ /Br  + (i  - f)th  g(n/Br] , (18) 

while  (ii)  efjr  has  to  obtained  from 

& = [/^/Br+(i-/)iin/Br]"1-[/p4ii/Br+ 

/(4nI-ifr)*^/Br  + (l-/)(^I-ifr)*^/Br]  • (19) 

These  two  equations  were  solved  iteratively  on  a computer.  The  obtained  numerical  results 
allowed  us  to  conclude  the  following  [2]: 

(i)  The  anisotropy  of  nonlinearity  can  be  considerably  different  from  the  anisotropy  of  the 
linearity  in  the  chosen  HCM. 

(ii)  Enhancement  of  nonlinearity  over  that  of  the  inclusion  material  is  possible,  the  enhance- 
ment being  anisotropic  too. 

In  closing,  we  note  that  the  algorithm  developed  can  be  easily  generalized  when  both  component 
materials  have  intensity-dependent  permittivity  scalars. 
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